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MICHAEL LUGO 

Abstract. This paper develops an analogy between the cycle 
structure of, on the one hand, random permutations with cycle 
lengths restricted to lie in an infinite set S with asymptotic density 
a and, on the other hand, permutations selected according to the 
Ewens distribution with parameter a. In particular we show that 
the asymptotic expected number of cycles of random permutations 
of [n] with all cycles even, with all cycles odd, and chosen from the 
Ewens distribution with parameter 1/2 are all | logro + O(l), and 
the variance is of the same order. Furthermore, we show that in 
permutations of [n] chosen from the Ewens distribution with pa- 
rameter a, the probability of a random element being in a cycle 
longer than approaches (1 — 7) " for large n. The same limit law 
holds for permutations with cycles carrying multiplicative weights 
with average a. We draw parallels between the Ewens distribution 
and the asymptotic-density case and explain why these parallels 
should exist using permutations drawn from weighted Boltzmann 
distributions. 



1. Introduction 

In this paper we study the cycle structure of random permutations 
in which the lengths of all cycles are constrained to lie in some infinite 
set S, and permutations may be made more or less likely to be chosen 
through multiplicative weights placed on their cycles. Cycle structures 
viewed in this manner are a special case of certain measures on S n 
which are conjugation-invariant and assign a weight to each element of 
S n based on its cycle structure. 

Definition 1.1. Let = (cr 1; <7 2 > • • •) be an infinite sequence of nonneg- 
ative real numbers. Then the weight of the permutation tt e S n , with 
respect to a, is 



w <?(tt) = Y\_ a 



Ci(-rr) 



i=l 



where c^tt) is the number of cycles of length i in tt. 
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Informally, each cycle in a permutation receives a weight depending 
on its length, and the weight of a permutation is the product of the 
weights of its cycles. The sequence a is called a weighting sequence. 

For each positive integer n, let (Q^,^) be a probability space 
defined as follows. Take = S n , the set of permutations of [n], 
and let F {n) be the set of all subsets of S n . Endow (fiW, F {n) ) with 
a probability measure P^ for each weighting sequence a as follows. 
Let P^(7r) = wJj(T)/^' e s ^Mi that is, each permutation has 
probability proportional to its weight. Extend P^ to all subsets of 
S n by additivity. To streamline the notation, we will sometimes write 
P^(tt) for F { j\ir). The sum of the weights of <j-weighted permutations 
of [n] is 

w <?{k) = n\[z n ] exp lj2°kZ k /k j 

7reS„ \fc>l / 

by the exponential formula for labelled combinatorial structures. 

We fix some notation. Define the random variable : — > Z + 
by setting X^ 1 \tt) equal to the number of /c-cycles in the permutation 

7i. Let X^(tt) = J2k=i -^l™^ 71 ") be the total number of cycles. We 
will often suppress 7r and (n) in the notation, and we will write (for 
example) F$(Xi = 1) as an abbreviation for F^({tt : x[ n \n) = 1}). 
Let Yfc = kXk- We define in order to simplify the statement of some 
results. 

This model incorporates various well-known classes of permutations, 
including generalized derangements (permutations in which a finite set 
of cycle lengths is prohibited), and the Ewens sampling formula from 
population genetics [8], which corresponds to the weighting sequence 
((7,(7,(7,...). If a is a 0-1 sequence with finitely many Is, then this 
model specializes to random permutations of which all cycle lengths lie 
in a finite set. These have a fascinating structure studied by Benaych- 
Georges [1] and Timashev [25] ; a typical permutation of [n] with cycle 
lengths in a finite set 5* has about j-n k / maxS fc-cycles, for each k in 
5*. In particular, most cycles are of length max S, which may be un- 
expected at first glance. Analytically, this situation is studied via the 
asymptotics of [^je^ 2 -* where P is a polynomial, as done by Wilf [28] . 
Yakymiv [29] has studied the case, alluded to by Bender [5] , in which a 
is a sequence of 0s and Is with a fixed density a of Is; the behavior of 
such permutations is in broad outline similar to that of the Ewens sam- 
pling formula with parameter a. An "enriched" version of the model 
has been studied by Ueltschi and coauthors [151 12S] . In their model, 



PROFILES OF PERMUTATIONS 



3 



permutations are endowed with a spatial structure. Each element of 
the ground set of the permutation is a point in the plane, and weights 
involve distances between points. Their "simple model of random per- 
mutations with cycle weight" [26j Sec. 2] is the model used here, where 
(Ti = e~ a \ 

There are other combinatorially interesting conjugation-invariant mea- 
sures on S n , including permutations with all cycle lengths distinct [T6j . 
and permutations with kth roots for some fixed k [TUJ [22]. However 
the generating functions counting these classes are not exponentials of 
"nice" functions and thus different techniques are required. 

Throughout this paper, we often implicitly assume that permuta- 
tions under the uniform measure on S n are the "primitive" structure, 
and weighted permutations are a perturbation of these. Here we follow 
Arratia et al. in [U, [2] , in embracing a similar philosophy and viewing 
the permutation as the archetype of a class of "logarithmic combina- 
torial structures", and Flajolet and Soria's definition of functions of 
logarithmic type [T%] . 

It will be convenient to use bivariate generating functions which 
count permutations by their size and number of cycles. In general, 
we take F(z,u) = ^ nk f n ,k^u k to be the bivariate generating func- 
tion, exponential in z and ordinary in u, of a combinatorial class J 7 , 
where f n ^ is the number of objects in T of size n and with a certain 
parameter equal to k. In our case n will be the number of elements 
of a permutation, and k the total number of cycles or the number of 
cycles of a specified size. Then [z n \ -j^F{z, u)\ _ 1 /[z n ]F(z, 1) gives the 
expected value of the parameter k for an object of size n selected uni- 
formly at random. The following lemma will frequently be useful, as it 
reduces the bivariate analysis to a univariate analysis. 

Lemma 1.2. Let f(z) be the exponential generating function of permu- 
tations with weight sequence a. Then the expected number of k-cycles 
in a permutation chosen according to the measure is 

Tff(ri) y Qk \z n - k ]f(z) 
9 k " k [*»]/(*) ■ 

Proof. The bivariate generating function counting the cycles of such 
permutations is 

o\z + 2 — H h Vk-i-j— J + u<J k-^ + Vk+i fc 1 H 

and this can be rewritten as (u — + Ylj>i ~~f~- Thus, from the 

exponential formula, the bivariate generating function counting such 
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permutations is 

V 3>l J J 

The expected number of cycles in a random permutation is [z n ]P u (z, l)/[z n ]P(z, 1), 
giving the result. □ 

The structure of this paper is as follows. In Section 2 we give exact 
formulas and asymptotic series (Propositions l2~2l and 12.31 ) for the mean 
and variance of the number of cycles of permutations chosen from the 
Ewens distribution. We also consider the average number of fc-cycles 
in such permutations of [n] for fixed k (Propositions 12.41 and 12.51) and 
for k = an (Proposition 12.61) . An "integrated" version of these results, 
Theorem 12 .7\ is one of the main results; this is a limit law for the 
probability that a random element of a weighted permutation is in a 
cycle within a certain prescribed range of lengths. In Section 3 we 
derive similar results for permutations in which all cycle lengths have 
the same parity. In addition, we determine the mean and variance of 
the number of cycles of such permutations (Theorem 13.61 treats the 
odd case, and Theorem 13.81 treats the even case). In Section 4 we 
explore connections to the generation of random objects by Boltzmann 
sampling. The main theorem of this section, Theorem I4.3[ states that 
the Boltzmann-sampled permutations of a certain class of approximate 
size n, including the Ewens and parity-constrained cases, have their 
number of cycles distributed with mean and variance approximately a 
constant multiple of logra. 



2. The Ewens sampling formula and Bernoulli 
decomposition 

The Ewens distribution [3J on permutations of [n] with parameter a 
gives to each permutation tt probability proportional to a x(jT \ This 
corresponds to the weighting sequence a = (a, a, a, . . .); we will write 
Pct , for Pl n ' ) ,El n ' ) , and call a random permutation selected in this 
manner a a -weighted permutation. In this section we derive formulas for 
the mean and variance of the number of cycles of permutations chosen 
from the Ewens distribution. Note that the number of cycles can be 
decomposed into a sum of independent Bernoulli random variables. 
Similar decompositions are due to Arratia et al. in J2j Sec. 5.2] for 
general a, and Feller 0, (46)] for a = 1; the fact that the number of 
cycles is normally distributed is seen in [HI Example 1]. Thus this 
section is largely expository; the proofs are provided for the purpose of 



PROFILES OF PERMUTATIONS 



5 



comparison with other proofs to be given below. The asymptotic series 
for and v£ n ^ appear to be new. 

Theorem 2.1. [20J Exercise 3.2.3] The distribution of the random vari- 
able X under the measure is that of the sum Y^k=i where the Zk 
are independent random variables and Z k has the Bernoulli distribution 
with mean o jio + k — 1 ) . 

Proof. The generating function of permutations of [n] counted by their 
number of cycles is Ylk=i S( n , = u{u + l)(u + 2) • • • (u + n — 1), 
where S(n, k) are the Stirling cycle numbers. Replacing u with o~u and 
normalizing gives the probability generating function for the number 
of cycles, 



and each factor is the probability generating function for a Bernoulli 



Combinatorially, we can envision this Bernoulli decomposition as fol- 
lows. We imagine forming a permutation of [n] by placing the elements 
1, . . . , n in cycles in turn. When the element k is inserted, with prob- 
ability ojio + k — 1) it is placed in a new cycle, and with probability 
1/ (cr+A; — 1) it is placed after any of 1, 2, . . . , k— 1 in the cycle containing 
that element. Then the probability of obtaining any permutation with 
c cycles is a c / (a(a + 1) • • ■ (a + n — 1)) , which is exactly the measure 
given to this permutation by P^. This is an instance of the Chinese 
Restaurant Process [201 Sec. 3.1]. 

From this decomposition into Bernoulli random variables, we can 
derive formulas for the mean and variance of the number of cycles 
under the measure In particular we note that since X is a sum 

of Bernoulli random variables with small mean, the variance of X is 
very close to its mean. Let ip denote the digamma function ip{z) = 
r'(z)/T(z)] this has an asymptotic series ijj{z) = log z — — j^z -2 + 
0(z~ 4 ) as z — >• oo. Let H n = Y12=l \ ^ e ^ ne harmonic number and 
let 7 = 0.57721 ... be the Euler-Mascheroni constant. 

Proposition 2.2. The expected number of cycles of a random a -weighted 
permutation of [n] is -X = a(tjj(n + a) — ip(cr)); in particular if a is 
a positive integer we have 



(1) E^X = alogn + (a 7 - aH a ^) + {a 2 - <7/2)n _1 + 0{n~ 2 ). 




random variable. 



□ 
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Proof. From Theorem 12.11 we have 



E ^x = y °- = aY 



i a + k-1 j-^a + k-l 

k=l fc=l 

Now, ip(z + 1) — ip(z) — 1/z; thus 

V»(n + a) - ij){a) = {^{n + a) - <0(n + a - 1)) H + {${o + 1) - ip{ 

1 1 1 
+ r + 



n + a — 1 n + cr — 2 a 



£ 

fe=i 



cr + fc - 1 



This proves that E„'X = a(^jj(n + a) — ip(cr)). The asymptotic series 
follows from that for ip(z) where we have used the fact that ip(n) = 
H n ^i — 7 when n is a positive integer. □ 

Proposition 2.3. The variance of the number of cycles of a random 
a -weighted permutation of [n] is 

(2) a 2 (V>'(n + a) - ij/(<r)) + a(ip{n + a) - if, (a)); 
this has an asymptotic series, 

A 2 1 

(3) V^ n) X = alogn + (-o-Y(o-) - <nf>(<r)) + °~ ~ rT 1 + 0{n- 2 ) 

The proof is similar to that of the previous proposition, noting that 
the variance of a Bernoulli random variable with mean p is p — p 2 . 

From ([3]) we can also derive for integer a the explicit formula (not 
involving ijj) 

cr+ra-l 

Y^X = -a 2 - + o-(\ogn + 1 -H a _ 1 ) + 0(l/n) 
which holds as n — > oo. It suffices to show that 

cr+n— 1 

(4) ^(n + ^)-^) = - E 

To see this, recall the identity ■?/>(£ + 1) — ^(sc) — 1/x; differentiating 
gives ip'(x + 1) — ip'(x) = —1/x 2 . Summation over x — a, a + l, . . . , a + 
n — 1 gives (jl]). 

Finally, we recall a normal distribution result for the total number of 
cycles [2, (5.22)]. Let X = x Jj£j?" be the standardization of X. Then 

limn^oo P^(X < x) — <&(x), where $(x) is the cumulative distribution 
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function of a standard normal random variable. This follows from 
Theorem 12.11 and the Lindeberg- Feller central limit theorem. 

We have thus far looked at the total number of cycles of a-weighted 
permutations. These distributions, suitably scaled, are continuous in 
the large-n limit. In contrast, looking at each cycle length separately, 
we approach a discrete distribution. More specifically, the number of 
&;-cycles of a- weighted permutations of [n], for large n, converges in 
distribution to V(a/k), where V(X) denotes a Poisson random variable 
with mean A; here we consider how quickly Eo -Xk approaches o~/k. 
Recall that Xk is a random variable, with Xk{it) the number of fc-cycles 
of a permutation tt. 

Proposition 2.4. [3], (37)]|27J The average number of k- cycles in a 
a-weighted permutation of [n] is 

(5) E^X k = ° } n)k 

k (n + a — l)k 

where (n)i~ = n(n — 1) . . . {n — k + 1) is the "falling power". 

We provide a new proof in terms of generating functions. 

Proof. The bivariate generating function counting cr-weighted permuta- 
tions by their size and number of /c-cycles is P(z, u) = (1— z)~ ua exp(er(u— 
l)z k /k). The mean number of /c-cycles is given by 

[z n ] d u P(z,u)\ u=1 _ [z n ]^(l - zy™ _ a [z n - k }{\ - 



[z n ]P(z, 1) [z n ](l-z)- a k [z n ](l-z)- a 

and the binomial formula gives (jSJ). □ 

When cr is an integer, a combinatorial proof can be obtained by 
considering a-weighted permutations as permutations where each cycle 
is colored in one of o colors. 

Proposition 2.5. There is an asymptotic series for E^Xk, 

Proof. The numerator and denominator of (jSJ) are polynomials in n 
of degree k; write the two highest- degree terms of each explicitly and 
divide. □ 

Proposition 2.6. Fix < a < 1. The expected number of elements in 
an-cycles of a random a-weighted permutation satisfies, as n —> oo, 

E^Y nn = ail - aY- 1 + 0( n - 1 ) 



s 
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(Here we have assumed for simplicity that an is an integer.) 

Proof. Let (3 = 1 — a. We have from Proposition 12.41 that 

, n) (n)an n\(Pn + a-l)\ n\ {(3n + a 

lhA >Y nn = a- — = a , - S1/ — ■ — T = o- 



{n + a-l) an ((3n)\(n + a - 1)\ (n + a - 1)\ ((3n) 

We now note that (n + r)!/n! = n r (l + 0(n -1 )), for constant r as 
n — > oo, from Stirling's formula. Applying this twice with r = a — 1 
gives the result. □ 

It would be of interest to determine the limiting distribution of the 
number of cycles with length between and 5n for constants 7 and 
5. There can be at most [7" 1 ] such cycles, so this random variable is 
supported on 0, 1, ... , |_7 -1 J • Thus to determine the limiting distribu- 
tion it suffices to determine the Oth through |_7 _1 Jth moments of this 
random variable. The a = 1 case will be treated in [T5] . 

We can essentially integrate the result of Proposition 12.51 to deter- 
mine the number of elements in cycles with normalized length in a 
specified interval. However, this can be done in a more general frame- 
work. Following [T3l Thm. VI. 1] and [T2J Thm. 1], for constants R > 1 
and <p > we define a A-domain as a set of the form 

A((f>,R) = {z:\z\<R,z=£ 1, |arg(>-l)| > <p}. 

Theorem 2.7. Let^2 k akZ k /k = a log j^- + K + o(l) be analytic in its 
intersection with some A-domain, for some constants a and K . Then 
the probability that a uniformly chosen random element of a random 
a -weighted permutation of [n] lies in a cycle of length between jn and 
5n approaches (1 — 7) " — (1 — 5) a as n — > 00. 

Note that analyticity in the slit plane suffices; this is the case = 0. 
We begin by stating two lemmas needed in the proof. 

Lemma 2.8. Let {a^'kLi be a sequence of nonnegative real numbers 
with mean a . Fix constants < 7 < S < 1. Then 

hm - T aJl--) = (l - 7 )" - (1 - 

fc=[ 7 n] v 7 

Proof. We rewrite the sum as an integral, 

L^ n J / 1 \ cr— 1 rSn / j \ o—l 



PROFILES OF PERMUTATIONS 



where fj,(x) = Ylj=i a j- Integrating by parts gives 

(6) (i - 5y- V(H - (1 - V(tw) - / K k )d ( 1 - - 



<5n / i \ cr-1 



70 



?? 



Differentiation allows us to rewrite the integral in ([H]) as a Riemann 
integral, 



(7) / ^(ifc)d 1 - - = - — - / dk) 1 - - dk. 



Let r(k) = fi(k) — ak. Then the integral on the right-hand side of ([7j) 
becomes 



a 1-- / r(fc) 1-- \ dk. 

n \./ 7n \ nJ J in V n ' ' 

We perform the first integral in (jSJ) and note that fi(5n) ~ cr-5n, [i{l n ) 

a ■ 7n in ([6]) . This gives 

(9) 

£ E f 1 - -V" - (i-7r-(i-*r+^ /%) f 1 - - x 

n k\ n ) n Jjn V " 

So it suffices to show that the final term in (Q is negligible, i. e. 



dk. 



/ T {k) 1 - - dA; = o(n 2 ). 
An V n J 



Since {cr^}^ has mean a, we have X}fe=i a k = cn + o(n). Thus r(/c) = 
o(n). On [771, 5n], (1 — k/n) a ~ 2 is bounded. So the integrand above is 
o(n), and the integral is o(n 2 ) as desired. □ 

Lemma 2.9. Say [z n ]P(z) = Cn cr_1 (l + o(l)) uniformly inn, for some 
positive constants C, a. Then 

^ fft [z-]P(z) ~ ^ H « 
as n 00, for any < 7 < 8 < 1. 

Proof. From the hypothesis that [;? n ]P(z) ~ Cn cr ^ 1 , we get 



[z n ]P(z) Cn*- 1 \ n 
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uniformly as n, k — > oo with < k < 5n. Therefore 

E-y^f - h (i+o(i) » 



5>H) +E— m-H) 



fc= |" 7 n] fc= |" 7 n] 

The first sum in the previous equation is O(n). The second sum has 
0(n) terms; since (1 — kjn) u ~ x and cr& can both be bounded above on 
the interval [771, Sn] each term is o(l). Thus the second sum is o(n). 
So 

x>y^ E(4-rn 

fc=r 7 n] fc=|" 7 n] \ v 7 / 

L5nj 



E H 1- s) (1+0(1)) 



= f 7 n] 

as desired. □ 



Proof of Theorem \2. 7[ This probability can be written as 

hm - > er fc - 

1— *nn n ' » 



Now, recall 



fc=|" 7 n] L J V 



^ a***/* = a log ^— + K + o(l) 



k 

by hypothesis. Thus the generating function P(z) of a- weighted per- 
mutations is 



'l-z)~ a e K (l+i 



P{z) = exp ^2 a kZ h /k^j = exp ^oTog + K + o(l)^ 

Applying the Flajolet-Odlyzko transfer theorem [T2], [z n ]P(,z) = Cn a ~ 1 (l+ 
o(l)) for some positive real constant C. Thus P(z) satisfies the hy- 
potheses of Lemma 12.91 Applying that lemma, we see that this sum 
is asymptotic to n~ x ^2k=Ln] a k(l ~ k/n) c ~ l \ the desired result then 
follows from Lemma [2.81 □ 

The hypotheses, and hence the conclusions, of Theorem 12.71 hold for 
many weight sequences o~x, 02, . . . with lim^oo - Ylk=i ° k = a > that * s ' 
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for weight sequences averaging a. In particular, we have the following 
special case. 

Corollary 2.10. Fix constants < 7 < 8 < 1. Let p a (n; 7 , S) be the 
probability that the element 1, in a a -weighted permutation of [n], lies 
in a cycle of length in the interval [yn, Sn] . Then 

lim Pa (n ]1 ,5) = (i- 1 y -(i-sy. 

Proof. We have the cycle generating function YlkLi azk /k — ^log tzi> 
apply Theorem 12. 71 □ 

For example, setting a = 1/2, 7 = 0.99, 5 = 1, we see that for large n, 
10% of elements of 1/2- weighted permutations are in cycles of length at 
least 0.99n. If we define the "co-length" of a cycle of a permutation to 
be the number of elements not in that cycle, a cleaner statement of the 
theorem becomes possible. The proportion of elements of a-weighted 
permutations in cycles of co-length at most (n is ( a . 

It would be desirable to replace the condition in the hypothesis of 
Theorem 12.71 with the less restrictive 

k 

it seems likely that this suffices to prove a limit law but the proof does 
not easily adapt to that case. 

3. Permutations with all cycle length of the same parity 

This section is devoted to results on random permutations in which 
all cycle lengths have the same parity; that is, they are either all even 
or all odd. We adopt the notation P^ for the family of measures 

where a = (0, 1, 0, 1, . . .), and similarly P^ for the family with a = 
(1,0,1,0,...); these are the measures corresponding to permutations 
with all cycle lengths even and with all cycle lengths odd, respectively. 

The results obtained here resemble those for the Ewens sampling 
formula with parameter 1/2. A heuristic explanation for this phenom- 
enon is as follows. Let us produce a permutation of [n] from the Ewens 
distribution with parameter 1/2 by first picking a permutation it uni- 
formly at random from S n , and then flipping a fair coin for each cycle 
of 7r. If all the coins come up heads we keep tt; otherwise we "throw 
back" the permutation 7r and repeat this process until we have a trial 
in which all coins come up heads. The number and normalized size 
of cycles of permutations obtained in this manner should be similar 
to those of permutations with all cycle lengths even, since for large 
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Figure 1. Ei 100) Y" fc for k = 2, 4, ... , 100 



permutations the parity constraint is essentially equivalent to a coin 
flip. 

Proposition 3.1. The expected number of elements in k-cycles of a 
permutation of [n] with all cycle lengths even is 

E (n) y n(n - 2) ■ ■ • (n - fc + 2) 

e k (n-l)(n-3)---(n-k + l) 

if k is even, and if k is odd. 

Proof. ByLemmaOl we have Y k = [z n - k }(l - z 2 )- 1 ' 2 /[z n }{\ - z 2 )- 1 ' 2 - 
we apply the binomial theorem and simplify. □ 

For example, when n = 10 we have 

(E( 10) F 2 , E( 10) F 4 , . . . , E( 10) F 10 ) = (10/9, 80/63, 32/21, 128/63, 256/63) 

w (1.11,1.27,1.52,2.03,4.06) 

and we observe that most entries are in the longer cycles. For n = 100 
this is illustrated in the figure above. 

Proposition 3.2. The expected number of elements in k-cycles of a 
random permutation of [n] with all cycle lengths odd is 



n(n—2) — (n-k+l) 
n-l)(n-3)-(ra-fc) ' 
ra-l)(n-3)-(n-fc+2) 
(n-2)(n-4)-(n-fc+l) ' 



n even 
n odd 
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Proof. The generating function of permutations with all cycle lengths 
odd, counted by their number of cycles, is P(z,u) = ((l+z)/(l — z)) u l 2 . 
We use Lemma [L2l to see that the mean number of elements in fc-cycles 

is given by E ( o ] Y k = [z n - k }^^/[z n }^^ z . We recall that [z r 

^ n ~^ V ' if n is even and n "^ n ~^ 2 " > " if n is odd; substituting and simplifying 
gives the result. □ 

By similar methods, we can obtain formulas for the exact number of 
permutations of [n] with all cycle lengths divisible by a, and the exact 
expected number of fc-cycles in such permutations for integers k which 
are divisible by a. These permutations have exponential generating 
function (1 — z a ) -1 / a . Permutations with all cycle lengths congruent to 
k mod a for some nonzero k are more difficult to deal with, as it appears 
that the generating function cannot be written in an elementary form 
except when a is even and k = a/2. (See [231 Sec. 5.0.3] for the relevant 
generating functions.) 

Proposition 3.3. (a) The number of elements ofk-cycles in a per- 
mutation of [n] with all cycle lengths odd, for fixed odd k, sat- 
isfies E^Yfc = 1 + + 0(n~ 2 ) as n approaches oo through 

even values, and Eo = 1 + 4^- + 0(n~ 2 ) as n approaches oo 
through odd values. 
(b) The number of elements of k- cycles in a permutation of [n] 
with all cycle lengths even, for fixed even k, satisfies Ee = 
1 + j- + 0(n~ 2 ) as n approaches infinity through even values. 

Proof. To prove (a), from Proposition 13.21 we have previous formulas 
for E^Yfc depending on the parity of n. These are fractions which have 
numerators and denominators which are polynomials in n; we can write 
out the two highest- degree terms of each polynomial and simplify. To 
prove (b) we proceed similarly from Proposition 13.11 □ 

Note that the expected number of elements in /c-cycles of permuta- 
tions with all cycle lengths even (or odd) approaches 1 as n gets large, 
if k has the appropriate parity. Assume we are dealing with permu- 
tations with all cycle lengths even. Naively, we might add the limits 
of the expected number of elements in 2-cycles, 4-cycles, . . . , n-cycles, 
and expect to get n. But these are each 1; their sum is n/2. Since each 
element is in a cycle, we must have J^LiE^Y^ = n. The difficulty is 

that the convergence ofE^n as n — ► oo is not uniform over k. Under 
the correct scaling, then, subtler phenomena can be seen; the "miss- 
ing" elements end up disproportionately in the longer cycle lengths for 
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permutations with all cycle lengths even. We note that similar phe- 
nomena of nonuniform convergence have previously been observed in 
random mappings, for example in [11] . 

Proposition 3.4. Fix e G (0,1). The expected number of elements 
in k-cycles in a random permutation of [n] with all cycle lengths even 
satisfies uniformly 

as k, n — > oo with < k/n < 1 — e. 

Proof. The result of Proposition 13.11 can be rewritten in terms of fac- 
torials as 



(n/2)! \ 2 (n-k)\ 



{{n-k)/2)\J n\ 
and by Stirling's approximation and routine simplifications, we have 



(10) E^Y k = (l - 



kY 1/2 l + i + 0(n- 2 ) 



Let n, k —>■ oo with < k/n < 1 — e. Then we have l/(4(n — 
k)) E [(An) -1 , (4en)~ 1 )], and so 0((n — A;)" 2 ) = 0(n~ 2 ). Furthermore 
l/(4(n — k)) = 0(n _1 ), with the constant implicit in the O-notation 
being (4e) _1 . Therefore 

/ £A~ 1/2 

E?»n= u--J (i + o^- 1 )) 

uniformly, as k, n — > oo with < k/n < 1 — e. □ 

Furthermore, we can essentially integrate the result of Proposition 
13.41 to determine the cumulative distribution function of the length of 
the cycle containing a random element of a random permutation with 
all cycle lengths even (or odd). This is the content of the next theorem. 

Theorem 3.5. Fix constants < 7 < S < 1. Let p e (n;^,5) be the 
probability that 1 is contained in a cycle of length between and Sn 
of a permutation chosen uniformly at random from all permutations of 
[n] with all cycle lengths even. Then 



lim p e (n; 7, 5) = v/I - 7 - Vl - 5. 

Since the measure P e is invariant under conjugation, this is the prob- 
ability that an element of [n] chosen uniformly at random is in a cycle 
of length between 771 and 8n in a random permutation of [n] with all 
cycle lengths even. 
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Proof. Note that 

£T = ^ log T^^ log rb + log2 + o(1) 

2\k 

and apply Theorem 12.71 □ 

The same is true for permutations with all cycle lengths odd; like 
those with all cycle lengths even they fall in the "cr = 1/2 class". 

We now move to consider the mean and variance of the total number 
of cycles of all lengths. 

Theorem 3.6. The mean number of cycles of a randomly chosen per- 
mutation of [n] with all cycle lengths odd is, as n — > oo, 

1 7 + 31og2 7 + logn /logn 
- log n H ± h O — ?r- 

2 & 2 8n V n 

where we take the + sign if n is odd and the — sign if n is even. The 
variance of the number of cycles is, as n — > oo, 

1 7 + 31og2-47r 2 ^log 2 n 
- log n H h (J 



2 ° 8 \ n 

Proof. We have the exponential generating function counting such per- 
mutations by size and number of cycles, (jz^) U ^ 2 . We can differentiate 
to obtain the mean and variance of the number of cycles. These are 
given by 



Let f r {z) = log r for r = 0, 1, 2 and let a r (n) = [z n }f r {z) for 

r = 0, 1, 2. Then we have 

a x (n) 2 a 2 (n) 2 

2ao{n) 4a (nJ 

and we need to find asymptotic series for the a r (n) as n — > oo. We 
observe that ao(n) is the number of permutations of [n] with all cycle 
lengths odd, which is (n — l)!! 2 /n! if n is even and n\\(n — 2)!!/n! if 
n is odd; Stirling's formula gives an asymptotic expansion, depending 
on the parity of n. To find a series for ai(n) as n — > oo, we expand 

fi( z ) = y fzf m a series with terms which are half-integral powers 
of 1 — z. From this we derive a series for J log y^f with terms of 
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the form (1 — z) 1 X I 2 V where L = log 1/(1 — z). The function being 
expanded is analytic in the complex plane slit along the real half-line 
{z e R : z > 1}; by [13, Thm. VI.3], an error of 0((1 -zf^^L) in the 
series for fi(z) leads to an error 0{n~ l ~ x l 2 logn) in the series for a\{n). 
We can thus transfer an asymptotic expansion for fi(z) near z = 1 
to give an expansion for ax(n) as n — > oo, and similarly for /^(z) and 
a 2 (n). Combining these series as specified by ffTTj) gives the result. □ 

We observe that this is (log 2) +o(l) more than the number of cycles 
of a permutation of [n] with all cycle lengths even. 

The following two results give a decomposition of the number of cy- 
cles of permutations with all cycle lengths even into a sum of Bernoulli 
random variables. 

Theorem 3.7. The generating function of permutations of [2n] with 
all cycle lengths even, counted by their number of cycles, is 

(12) p 2n (u) = [u(u + 2){u + 4) • • - (u + (2n - 2))] ■ (2n - 1)!! 

Proof. The bivariate generating function for permutations with all cycle 
lengths even, counted by their size and number of cycles, is (1 — z 2 )~ u l 2 . 
Let pk(u) be the desired generating function. Then we have 

z 2 

(1 - z 2 y u/2 = po(u) + pi{u)z + p 2 {u)— + ■■■ 

and it is clear that pk is the zero polynomial for odd k. For even k, the 
binomial theorem gives 

(i - z 2 r /2 = i + (~ u / 2 ) (-^ 2 ) + (~ U J 2 ) i-z 2 f + ■■■ 



1 j \ 2 

and so we have P2n(u) = (2^)'( by comparing coefficients; this can 
be expanded to give the expression above. □ 

A combinatorial proof is also possible. Recall that we can write a 
permutation 7r of [n] in terms of its inversion table, a sequence of inte- 
gers a%, ci2, • • • , a n , with a, = \{j : j < i, ir(j) > 7r (z) } | . The number of 
zeros in the sequence (a 1; . . . , a n ) is the number of left-to-right maxima 
of 7T. The "fundamental correspondence" between permutations writ- 
ten in cycle notation and in one-line notation takes permutations with 
k cycles to those with k left-to-right maxima; furthermore, permuta- 
tions with all cycle lengths even are taken to those with all left-to-right 
maxima in odd positions, and conversely. Thus it suffices to show that 
p n (u) is the generating function of permutations with all left-to- right 
maxima in odd positions, counted by their number of maxima; this is 
done by considering the inversion table. 
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Theorem 3.8. The number of cycles C n of a random permutation of 
[2n] with all cycle lengths even, as n — ► oo, is asymptotically normally 
distributed with 

(13) E(C n ) = Uogn+ Q 7 + l 0g 2j +0(n- 1 ) 
and 

(14) V(C n ) = l - \ogn + Q 7 + log2 - ^] + Oin- 1 ) 

Proof Let n = 2m. From Theorem 13.71 we have C m = Y^k=i X m ,k 
where the X m ^ are independent Bernoulli random variables with F(X m ^ 
1) = 1/ (2k — 1). The formula (1131) for the expectation follows from the 
asymptotic series for the harmonic numbers. The variance is given by 

m / /1\ 2 \ - 1 

and we need to consider the second sum. We have Y^rj=\ = —^'( m + 
1) + tt 2 /6, so 

/ 1 7T 2 ' 

YC m = EC m - f -V>'(2n + 1) + -if/{n + 1) + — 

But ^'(m) = O^" 1 ), so in fact we get YC m = EC m -f + 0{l/m); 
from this and f|T3|) we get f|14l) . Asymptotic normality follows from the 
Lindeberg- Feller central limit theorem [7]. □ 

There is not such a simple decomposition for permutations with all 
cycle lengths odd. However, it appears that the polynomials counting 
permutations of [n] with all cycle lengths odd by their number of cycles 
have only pure imaginary roots. If this is true, then the number of 
cycles of a random permutation of [n] with all cycle lengths odd can 
be decomposed into a sum of [n/2\ independent {0, 2}- valued random 
variables, plus 1 if n is odd. It may be of interest to study the zeros of 
these polynomials. 

4. BOLTZMANN SAMPLING 

We have at this point seen substantial similarities between permuta- 
tions with all cycle lengths having the same parity and permutations 
with cycle weights 1/2. This suggests that an average of weights is in 
some sense a more fundamental parameter than the individual weights. 
This has been anticipated by the notion of a function of logarithmic 
type [H], which has been used in the study of permutations \17\. Let 
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A (p, rj) = {z : \z\ < p + rj, z £ [p, p + rf\}. A function G(z) is called 
logarithmic if it is of the form 



for some constant multiplier a and function R(z), where R(z) is ana- 
lytic in A and satisfies R(z) = K + o(l) for some constant K as z — > p 
in Ao, and p is the unique dominant singularity of G on its circle of 
convergence. In [H, Prop. 1] structures having components enumer- 
ated by a function of logarithmic type G(z) are considered; for such 
structures of size n, the expected number of cycles is alogn + 0(1), as 
is the variance. However, the structures considered in this paper have 
not all had components counted by functions of logarithmic type. For 
example, the components of permutations with all cycle lengths even 
are counted by the exponential generating function ^log^a, which 
has singularities at z — ±1 and thus does not have a unique dominant 
singularity. 

The following conjecture, in the light of these averaging phenomena, 
seems natural. It is supported by Theorem 14. 3[ an analogous result on 
"Boltzmannized" permutations. 

Conjecture 4.1. Let a = (o~i, o"2, . . .) be a sequence of nonnegative real 
numbers with mean a, that is, with lim^oo - (%2k=i a k) = a - Then 
permutations of [n] selected according to the weights a have an asymp- 
totically Gaussian number of cycles asn — >• oo, with mean and variance 
asymptotic to alogn. 

Definition 4.2. Let a = (01,02, . . .) be a weighting sequence, and let 
x be a positive real parameter. Let \ir\ denote the size of the ground set 
of a permutation tt. Then we define the a -weighted Boltzmann measure 
with parameter x on permutations, a probability measure on [J^ =0 Sk, 
by 



(See Definition 11.11 for the weight wg(it).) 

For any choice of a and x, Fff tX is a probability measure. It suffices to 
show that has total mass 1. But J2k>i <J k xk /k is the weighted gen- 
erating function of cycles, and we can apply the exponential formula. 
Thus ¥ff )X is a straightforward weighted generalization of the Boltz- 
mann measure on labelled objects studied in [HJ [10]. We also retain the 
formulas from [6j Thm 2.1] for the expected size and the variance of 
the size of the objects chosen according to this measure. Let C(x) be 






ex P (£*>! °kX k 



/k) 
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the exponential generating function of a labelled combinatorial class, 
and N the size of a random object chosen from that class according to 
the Boltzmann measure. Then 

E- (N) = X ^ C{ ~ X) E- (N 2 ) = ^ 2 °^ 
aA V C{x) ' a ' x[ ' C(x) ■ 

We now assemble a sequence of lemmas. These lemmas will be used 
to prove the following theorem, which is the main result of this section. 

Theorem 4.3. Let a = (o~i, 02, . . .) be a weighting sequence with mean 
a. Let x = x(fi) be chosen so that K^ yX (N) = fi. Let X be a ran- 
dom variable denoting the number of cycles of a permutation. Then 
Eff >x (X) = Yff >a: (X) ~ cdog/i as x — > l~ or ji — > oo. 

The main analytic result needed follows. 

Lemma 4.4. [2TJ Exercise 1.88] Let bo, b±, . . . be positive real numbers, 
such that Yl^Lo °n is divergent, and J2k>o ° k ^ k is convergent for < 
t < 1. Then 

a + aH ha„ . Y,k>o a kt k 

hm — = s implies hm = — - = s. 

bo + oi H Vo n t-+i- 2^ fc>0 M fc 



Lemma 4.5. Let cj., 02, (73, . . . be a sequence of real numbers such that 

n 

, a k = a. 



1 - 
lim — \ 1 

n— >oo 77, ^— — ' 

k=l 

Then Er=i **** = + - x)- 1 ). 

Proof. Apply Lemma [4.41 with a k = o~k,bk = 1. □ 

Lemma 4.6. Let {o~k} k % 1 be a sequence of nonnegative real numbers, 
bounded above, such that Ylk=i ak ~ an as n ~ ¥ 00 > f or some constant 
a > 0. T/ien Ylk=i it ~ a ^°S n as n ~ ¥ 00 ■ 

Proof. We begin by showing that if J"™ /(x) rfx ~ n as n — > 00 for some 

function / such that f{x)/x is integrable on [1, 00), then J™ dx ~ 

logra as n — > 00. Let = f(x) dx. We integrate J™ ^— dx by 

parts, getting 



x n 1 J 1 x 2 

Clearly F(l) = 0, and F(n) ~ n by assumption, so 



Jl x J 1 X 2 
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Since F(x) ~ x as x — > oo, the integrand satisfies F(x)/x 2 ~ 1/x, and 
so 



dx 



ar 



n j 



<ix = log n, 



x 



proving the claim. 

Now, we need to check that this statement about integrals translates 
into an analogous one about sums. Let {ak}^ =1 be as in the hypothesis, 
and let f(x) = a^ x \. Then we want to show that 



n+l 



/(*) 



X 



dx-J2 



o(logn) 



as n —>■ oo. We have 



n+l 



/(*) 



X 



dx — 



/(*) 



fe=i 



fe=i 



E 

fc=i 



fe+i 



/(*) 



,r 



lM log HB) 



and so, since | log(l + l/k) — l/k\< I /2k 2 for positive integer k, 



n+l 



/(*) 



dx — 



/(*) 



fe=i 



k 



< 



^ 1k 2 



k=l 



Since {f(k)}^ =1 is bounded, the sum on the right-hand side is conver- 
gent. We have f™ +1 f(x)/xdx ~ logn, so Ylk=i /WA ~ logn as well. 
Thus we have proven the lemma for a = 1. Multiplying through by a 
gives the desired result. □ 



Lemma 4.7. Let o~i, a% 03, . . 
such 

lim^oo J ]T™ =1 = r ^ ere 



fee a sequence of positive real numbers 

that 



(15) 



OO £, 



E 

fc=i 



a log 



1 — x 



o log 



.r 



Proof. Applying Lemma 14.61 to the hypothesis, lim^oc j^L_ Y^™ =1 2k 
a. We apply Lemma S3] with a k = o~k/k,bk = l/k. This gives us 

lim EL^A = lim Sjsfffi. 

n-»oc l+if n x->l- 2^ k>1 X k /k 
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Now, ^2 k>1 x k /k = 1°S(V(1 — x ))> an d 1 + ^» ~ logn, so we have 

J2k>l a kX k /k 



lim N — = lim 

ra-»oologn'^ A; log(l/(l — x)) 

Thus the right-hand side here has value a, proving (|T5|) . □ 

Proof of Theorem \4-3[ Note that for the Boltzmann measure with pa- 
rameter x and weight sequence a, we have C(x) = exp w2 k>1 o~kX k /k). 
Thus the distribution of sizes N under this measure has expectation 

\fc>i / fc>i 

Furthermore, the Boltzmann distribution P^ can be obtained by tak- 
ing V(o~kX k /k) cycles of length k for each k > 1 and forming uniformly 
at random a permutation with the resulting cycle type. The mean and 
variance of the number of cycles chosen from the distribution P^ is 
thus exactly J2k>i °~kX k /k. Since we have E-^iV) ~ a/ (1—x) as x — ► 1~ 
by Lemma H31 we can solve for x to see that 1 — E ~ x as x — ► 1~. 
Therefore 

o~kX k / k ~ a log ~ a log 7 r- = a log — — — - ~ a log E x (iV) 

which is the desired result. □ 



It would be desirable to translate Theorem 14.31 into a result about 
permutations of a fixed size selected uniformly at random; this is one 
possible way of proving Conjecture 14.11 Note that P^ iX is a mixture of 

the various . It is often possible to prove results about a family of 
measures P\, parametrized by A, which are mixtures of well-understood 
measures P^, where we draw from P*™** with probability e~ x X n /n\; 
this goes by the name of analytic de-Poissonization [18], [24"l Ch. 10]. 
Informally, we pick from P^ where N is Poisson with parameter A. 
In the case described here we can get results on permutations chosen 
from P^) where N is the size of objects from a Boltzmann distribution; 
thus techniques of "analytic de-Boltzmannization" will be necessary to 
achieve this goal. 
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